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A B-spline  Tensor  for  Vectorial 
Quasi-Interpolant 


Fabrice  Dodu 


Abstract.  The  aim  of  this  paper  is  to  introduce  new  techniques  and 
new  tools  for  vector  field  approximation.  We  do  so  by  building  the  equiva- 
lent of  B-splines,  which  are  now  tensor  B-splines,  as  shown  below,  and  by 
appling  to  the  discretization  based  on  a regular  grid  of  a differential  op- 
erator a fundamental  solution  of  it,  as  done  for  polynomial  B-splines  and 
polyharmonic  B-splines  (see  [5,6]).  We  thus  obtain  quasi-interpolants  in 
the  vectorial  case  whose  properties  generalize  the  properties  of  the  quasi- 
interpolants  generated  by  using  B-splines.  All  this  is  done  in  the  case 
when  the  data  lie  on  a regular  infinite  grid. 


§1.  Introduction 

Fluid  mechanics,  meteorology  and  more  and  more  other  applications  need 
approximate  functions  from  1R3  to  R3.  Given  some  discrete  vectorial  data, 
we  want  to  get  a function  interpolating  or  approximating  the  data.  At  first 
glance,  we  may  think  of  doing  this  with  three  independent  approximations 
of  the  data  (one  for  each  component  of  the  data).  Of  course  this  can  be 
done,  but  it  usually  gives  poor  results  since  there  is  no  connection  between 
the  various  components  of  the  approximation  function,  while  the  applications 
may  require,  for  example,  a divergence-free  (or  a rotational-free)  function.  In 
order  to  take  into  account  this  kind  of  connection,  we  want  to  determine  the 
function  interpolating  the  data  and  minimizing  a seminorm  over  all  vectorial 
functions  interpolating  the  data.  The  seminorm  which  is  minimized  is  based  on 
the  Helmholtz  decomposition  of  vector  fields  into  a rotational  and  a gradient 
part  ( p ||div  • ||2  + 1 1 rot  • |j2).  This  will  be  presented  in  detail  in  the  forthcoming 
thesis  [2],  The  weight  p is  introduced  to  allow  the  rotational  part  of  the  field 
to  dominate  the  gradient  part  and  so  in  the  so  obtained  function  [2], 

In  order  to  determine  the  interpolating  vector,  we  need  to  solve  a linear 
system  which  is  usually  large  (three  times  the  number  of  data)  and  badly 
conditioned.  This  is  why  in  this  paper  we  do  not  propose  to  interpolate, 
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but  instead  approximate  the  data  by  building  a quasi-interpolant  based  on 
5-splines.  Note  that  if  (zj)"=1  are  vectors  and  if  we  want  to  get  a vectorial 
function  S such  that  S{x)  — 5;(x)  z,-  (as  we  do  in  the  scalar  case  with  5- 

splines  5,),  we  need  that  the  functions  Bi  are  matrix  and  not  scalar  functions 
as  in  the  case  of  scalar  approximation  in  It  or  Rd.  This  is  why  the  tools  we 
build  are  tensors. 

Now,  in  order  to  build  this  “5-spline  Tensor”,  we  will  use  the  same  strat- 
egy as  for  polynomial  and  polyharmonic  5-splines:  we  will  discretize  the 
differential  operator  P7n,p  (5)  defined  in  [2]  and  apply  to  this  discretization  a 
fundamental  solution  of  Pm,p  (5) , thus  obtaining  a kind  of  approximation  of 
the  Dirac  tensor  6 1 3.  All  this  is  done  on  a cardinal  grid  (i.e.  a regular  infinite 
grid). 

For  our  three-dimensional  problem  [2],  we  choose  the  3x3  differential 
matrix  defined  by 

Pm,p  (5)  = (-l)m  Am_1  [pVdiv  • —rot  (rot-)] , 

where  div  is  the  divergence  operator,  V is  the  gradient  operator,  rot  is  the 
rotational  operator  and  p is  an  arbitrary  positive  parameter. 

We  remark  that  if  p = 1,  then  Pm  1 ( D ) = (— l)m  A 771  J3,  where  1 3 is  the 
identity  matrix  of  R3,  so  we  obtain  three  independent  operators  (one  for  each 
component),  each  one  being  as  in  ([5]). 

In  the  first  part  of  this  paper  we  give  the  construction  of  a discretization 
of  Pm,p  (5).  In  the  second  part,  we  define  polyharmonic  5-spline  tensors  and 
give  the  main  properties.  In  the  last  part,  we  study  the  associated  vector 
quasi-interpolants,  and  in  particular  prove  their  P^-reproduction. 


Notation:  Let  m be  a integer  with  m>  2.  P„  denotes  the  set  of  polynomial 
with  variable  in  R3  of  total  degree  at  most  n and  Vn  = Pn  x Pn  x Pn.  P> 
denotes  the  set  of  distribution  of  R3  and  T>  = D x D x D . Let  / : R3  1 — + R 
be  a scalar  function  of  a three  dimensional  variable.  Let 


( £>2/3  - £>3/2  \ 
rot/  = d3fi  -di  f3 
\d1f2  — £>2/1  / 


be  the  rotational  operator,  and  let  div/  = £>1/1+02/2  + 03/3  be  the  divergence 


operator.  Let  A771  = 


( x^3 


.*)’ 


T (g)  is  the  Fourier  transform  of  the 

function  (or  distribution)  g.  For  all  ( € R3,  sin  (£)  denotes  the  vector  defined 
by  VI  < i < 3,  sin(£)j  = sin(£j).  6 denotes  the  Dirac  distribution.  Let 

be  the  function  such  that  : vm+i  = ^1+7^3/ 2 II  ■ ||2m_1-  All  tensors 
will  be  denoted  with  bold  capital  letters  (i.e.  X,...).  Let  h > 0,  and  1 be 
such  that  : i2  = — 1.  (e,)3=1  denotes  the  canonical  basis  of  R3.  For  k € IN 
and  i = 1,2,  or  3,  6^  denotes  the  fcth  divided  difference  of  step  h defined 


s{h,i  1]  0 6h,i ■ We  use 


by  (6h,i /)  (*)  = / (*  + te<)  - / (* - le0-  6h,i 

standard  multi-index  notations.  If  a G IN3  and  |a|  = Oj  + Q2  + 0:3,  then 
Da  = 


jri 


8°2 


d° 3 


(9xi)“!  [dx2)a2  (Sx3)“3  ’ 


„Ol„«2_«3 

L1  x2  x3  ■ 
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§2.  Discretization  of  Pm<p(D) 

The  goal  of  this  section  is  to  define  an  approximation  of  Pm,p  (D)  which  should 
reproduce  polynomials  of  largest  possible  degree  (in  the  sense  defined  below). 
We  now  introduce  the  P^-exactness  of  an  operator,  and  a Pfc-exact  approxi- 
mation of  Da. 

Definition  1.  An  operator  E'  is  said  to  be  a P^-exact  (resp.  Vk-exact)  ap- 
proximation of  an  operator  E iff  for  any  function  f ,E'  f is  a linear  combination 
of  translates  of  f and  for  any  p € P*,  (resp.  p (zVk),  E' p = Ep. 

Remark.  h~1d\  i is  a P 2 -exact  approximation  of  jP. 

Definition  2 . Let  D\  < N be  the  operator  defined  by 


D1 


N 

E 

fc= 0 


(-1)*  (*02  (cl  nc2M 

~~2h~  (2fc"+T)!  i 


For  any  a € IN3,  we  define  an  approximation  of  step  h and  level  N of  Da  to  be 
the  operator  defined  by 


D°h>N = (DiXNr  o (DU,Nr  ° (du 


3 

n)  ■ 


In  the  same  way,  Pm,p  ( Dh,N ) denotes  the  approximation  of  Pm,p  (D). 

Proposition  3.  Let  a € IN3  and  D%  N be  defined  as  above.  Then,  for  all 
mappings  f from  1R3  to  IR,  there  exists  real  constants  ( c7 ) such  that 

F (Dh,Nf)  = ( i E c7  (sin  W)7 ) r (/)  • 

V |7|=(2AT+1)M  / 

Proof:  For  every  1 < i < 3, 

? ( sh,if ) (C)  = (exp  (i7r/i£j)  — exp(—inh<(i))  ■ IF (f)  (C)  = 2isin(h()i  ■ E (f)  (C) 

By  applying  the  Fourier’s  transform  to  the  (2 N + 1)  |a|  centered  differences 
of  D'l  N,  we  obtain  the  result.  □ 

Proposition  4.  Let  a € IN3.  Then 

i)  _D£  N is  a P|a|+2JV+i'exact  approximation  of  Da . 

ii)  Pm,p  ( Dh,N ) is  an  approximation  of  Pm,P  (D)  which  is  V2m+2N+i-exact. 

Proof:  We  prove  by  induction  on  |a|  that  for  all  / € C2W+2+l“l , there  exist 
c7  real  constants  such  that 


DZNf  = D“f+  J2  c7  DViCy). 

\-y\=2N+2+\a\ 
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In  the  following,  for  any  7 £ IN3,  dy,  e7  are  real  constants.  The  proof  for 
|a|  = 1 is  due  to  Steffensen  in  ([7]).  Suppose  that  it  is  true  for  |a|  = m.  Let 
|a|  = m + 1,  / £ (72JV+2+m+i  arKj  jej.  us  c}j00se  j such  that  a,  0.  Then 

Dah,Nf  = DlitN  (Dah^f) 

= DlitNlDa-^f+  £ cyDV{  Cy) 

\ |7|=27V+2+m 

= D}D°-*‘f+  Y1  dyDT{Da-'-f(Cy))+  Y,  c^7/(C 7) 

|7|=2W+3  |7|=2iV+3+m 

= Daf+  Y e-r/(C7)- 

|'y|=2iV-f-3-{-m 

Thus,  we  get  the  result  if  we  note  that  Pm,p  ( D ) is  a differential  operator  of 
degree  2m.  □ 

§3.  B-spline  Tensors  Associated  with  Pm  p(D) 

Lemma  5.  Let  vm+\  be  such  that  (— l)m+1  Am+1vm+i  = S in  D' . Let 

x = tDl  v2um+1 

p 

({p2,2  "I"  ^3,3)  1 — 9j  2vm  + l — @l,3vm+l  \ 

— 01,2»m+l  (<^,1  + ^3,3)  ^m+1  -0|i3Vm  + i J . 

-^l,3?m+l  -5|j3Vm+i  (dh  + dl  2)  Vm+l  / 

Then  PmiP{D)X  = 6I3. 


Definition  6.  Let  X be  a fundamental  solution  tensor  of  Pm,p  ( D)-X  = 5-I3. 
We  define  the  level  N and  step  h S-spline  tensor  associated  with  the  operator 
Pm,p  ( D ),  to  be  the  tensor  B™N  p defined  by  = h3  PmiP  (Dh:N)  X. 

Remarks. 

a)  B™’£  is  not  a symmetrical  tensor. 

b)  Ifm  is  even,  we  can  prove  the  existence  of  a differential  matrix  Rm.,p  (D) 
of  degree  m such  that  : Pm,p  (D)  = R™  p{D)  RmtP{D)  and  we  can 
construct  a symmetrical  tensor  C^N  p defined  by 


C 


m,p 

h,N 


= h3  R^,P  {DhtN)  X Rf,p  {Dh,N)  . 


Lemma  7.  The  elements  of 
tributions  on  1R3. 


l<tj,<3 


are  in  the  set  of  tempered  dis- 
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Denoting  by  T the  Fourier  transform  applied  to  each  element  of 

the  tensor  , we  obtain  the  following  theorem: 

Theorem  8.  With  the  above  notation, 

i)  (0)  = h3-I3, 

ii)  V7  € IN* , 111  < 2N  + 1 ; D1T  (b?$  (0))  = 0, 

iii)  VI  <j,l<  3,  te).  (t)  = 0(\t\-™~5), 

V 'Ji*  |<H+oo 

i»)  VJs  £ 7Z?„  VT  £ WJ,  |7|  < 2m  - 1 ; WT  (b~,^  (f)  = 0, 

v)  ^5  = B5f(i). 

Proof:  For  i)— ii),  we  use  Definition  2 and  Proposition  3.  By  the  Taylor 
expansion  of  sin,  near  0 we  obtain  T (-8™$)  (0  = h3  I^  + Q (£),  where  Q (Q 
is  such  that  for  all  i and  j such  that  1 < i,j  < 3,  there  exist  real  coefficients 
such  that  Qt  J (C)  = E|a|>2N+2m+4CaJ'pi^pTFT-  Thus>  we  obtain  the 
results. 

We  prove  iii).  By  using  the  above  expression  of  Q,  D1  T are 

integrable  near  0 if  and  only  if  2 N + 2 — |i|  > —3.  Using  Proposition  3 
and  Lemma  5,  there  exist  real  coefficients  d'^’k  and  for  all  i and  j such  that 
1 < i,j  < 3 


(KS),/ o = E 

J k= 1 


1 

h <?m 


E 


|a|=2m(2JV+l) 


(sin(7rhC))a  t f 
||27rC||2m+2 


As  a consequence,  the  elements  of  T [B™^ j 
nal  fractions  of  degree  —2m.  D^T 
only  if  -2m  - |7|  < -3.  Then,  D7f 


are  bounded  at  infinity  by  ratio- 

are  integrable  near  infinity  if  and 

are  integrable  in  It3  iff  2N  + 
i.i 


2 - |7|  > -3  and -2m- |7|  < -3.  Thus,  Wp(b?£)  = o (\t\-2N~4), 

' ’ /j,i|t|-»+oo 

for  all  1 < j,  l < 3.  Using  [4],  the  last  expression  may  be  strengthened  to 


ZUJF 


= 0(|t| 

|t|-+oo  vl 


-2JV-5 


). 


1 <j,l  < 3. 


We  now  establish  iv).  For  all  j in  72? , j in  IN3  and  positive  h,  we  have 
D1  T (i)  = 0 iff  D~*T ( Pm,p  {Dh,N))  (0)  = 0.  Now  by  using  Proposi- 

tion 3,  we  have,  for  all  tensors  Y and  all  1 in  IN3  such  that  |7|  < 2 N + 2m+ 1: 


d7  (r{pm,p  (Dh}N)  y))  (0)  = d7  [r{pm,p  (D) ) (o)  ■ t (y)  (o)  ) . 
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Furthermore,  T ( Pm,p  ( D ))  is  a polynomial  matrix  of  degree  2 m so  V|-y|  < 
2 m — 1 , D7F  ( Pm,p  ( D ))  (0)  = 0.  This  gives  iv). 

Finally,  we  prove  v).  According  to  the  Fourier  transform’s  properties,  v) 
is  equivalent  to 

r {Kn)  (0  = T {Kn  (x))  (0  - h3  ■ r [Kn)  (h  ■ 0 . VC  e ]R3. 

Now,  from  the  definition  of  B™$,  we  derive  for  any  C in  H3, 

T (b£#)  (ft  • C)  = ^ (Pm,p  (BliJV)  *)  (ft  • C) 

= T(pm,p  ) (ft  • C)  • T (■ Pm,p ) (ft  • C)'1  • 

Using  Proposition  3,  we  have 

T{Pm,P  (D1}N))  (ft  • C)  = h2mT(Pm,p  {Dh,N))  (C) , 

F(Pm,P  (D))-1  (ft  • 0 = h-2mf(Pm,p  (D))"1  (C) , 
and  thus  VC  € 1R3, 

P (*££)  (h  ■ C)  = P(Pm,p  (D*,jv))  (C)  • P(Pm,p  (^))_1  (0 

= ft-3-P(B^)(C).  □ 


Remarks. 

a)  These  polyharmonic  B-spline  tensors  may  be  considered  as  a regularisa- 
tion  of  the  Dirac  distribution  tensor  61 3. 

b)  We  obtain  the  same  properties  with 


§4.  Associated  Vector  Quasi-Interpolant 

Given  vectorial  data  (uOjezz3,  in  this  section  we  define  a vector  field  S ap- 
proximating the  data  (jh,  Zj)^e7/p  (he.  such  that  S(jh)  ~ Zj  for  all  j in 
77j  ‘‘),  by  using  the  above  defined  tensor  B-splines.  This  vector  generalizes  the 
polyharmonic  B-spline  quasi-interpolant  (see  [5,6]). 

Definition  9.  Let  B™'£  be  the  level  N and  step  ft  B-spline  tensor  associated 
with  P-m.p  (D).  For  all  j € 7Z3,  let  Zj  (E  R3,  and  let  z = (zj)^e2Z 3-  Then 
the  vector  quasi-interpolant  of  step  ft  and  level  N associated  with  the  operator 
Pm,p  (D)  and  the  {jh,Zj)-e7L 3 data,  is  the  vector  function  defined  by 

zp 

1? 
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Theorem  10.  Let  l = inf  {2Ar  + 1, 2 m — 1},  and  suppose  there  exists  p £Vi 
such  that  for  all  j in  ZZ3,  Zj  = p{jh).  Then,  = p.  We  say  that  the 

vector  quasi-interpolant  of  step  h and  level  N reproduces  Vi-  As  a particular 
case,  S™£_1  reproduces  V^m-i- 

Proof:  The  proof  is  based  mainly  on  Poisson’s  egality  and  Theorem  8.  It 
follows  along  the  same  lines  as  the  proof  in  [3]  in  the  scalar  case.  □ 

Remarks. 

a)  2 m — 1 is  the  maximal  order  of  reproduction. 

b)  We  obtain  the  same  properties  if  we  define  the  vector  quasi-interpolant 
using  by  the  symmetrical  tensor  C™'f, . 

c)  A similar  problem  is  studied  in  H2  in  ([1]),  where  using  another  discretiza- 
tion of  P2,p(D),  the  authors  obtain  a vector  quasi-interpolant  which  is 
(Pi  (IR2))  -reproducing. 

Theorem  11.  Let  f be  a vector  function  of  C (lR3)k -class  and  all  partial 
derivatives  of  order  k being  bounded  over  IR3.  Let  S™$.j  be  the  above  defined 

vector  quasi-interpolant  associated  to  the  (jh,  f ( jh ) ^ data.  Then 


sup  || S™$.f  (t)-/(t)||  = 


t€  R3 


h—0 


0(hk) 

0(h2N+2\ln{h)\) 


ifk  < 2JV  + 1, 
ifk  >2N  + 2. 


Proof:  The  proof  follows  that  of  Theorems  4.11,  5.1  and  5.6  in  [3].  □ 
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